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Abstract 



The inclusive production of charged hadrons in the collisions of quasi-real photons 
(e + e~ — > e + e~+X) has been measured using the OPAL detector at LEP. The data were 
taken at e + e~ centre-of-mass energies from 183 to 209 GeV. The differential cross-sections 
as a function of the transverse momentum and the pseudorapidity of the hadrons are com- 
pared to theoretical calculations of up to next-to-leading order (NLO) in the strong cou- 
pling constant a s . The data are also compared to a measurement by the L3 Collaboration, 
in which a large deviation from the NLO predictions is observed. 
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1 Introduction 



Hadronic interactions of two photons lead to the production of hadrons whose properties in 
Quantum Chromodynamics (QCD) depend on the underlying partonic processes, as well as on 
the way in which the partons are transformed into observable hadrons. The study of inclusive 
hadron production is therefore an appropriate tool to investigate the validity of QCD in hadronic 
photon-photon interactions. Furthermore a drastic discrepancy between data and NLO QCD 
has recently been observed in this process [1] which would indicate the breakdown of this 
otherwise so successful theory and hence requires further study. Our measurements complement 
similar studies of jet production [2,3]. The same partonic processes are at work in QCD in 
both cases, but hadronisation is usually treated differently. To predict the production of single 
hadrons, universal fragmentation functions obtained from fits to independent data are utilised to 
describe the transition from partons to hadrons. Jet observables on the other hand are designed 
to keep the influence of hadronisation processes small. The remaining distortion with respect 
to partonic jet observables is estimated using hadronisation models. The reliable calculation of 
an observable in perturbative QCD requires a sufficiently high energy scale which in our case is 
provided by the transverse momenta of the charged particles. First results have been published 
by OPAL at e + e~ centre-of-mass energies y/s^ = 161 and 172 GeV [4]. L3 has published 
results in the range y/s^ = 189 — 202 GeV [1]. This paper extends the OPAL measurements of 
the transverse momentum spectra and pseudorapiditjo distributions of charged hadrons in [4] 
using data taken at e + e" centre-of-mass energies from yfs^, = 183 to 209 GeV, representing a 
roughly thirty-fold increase in integrated luminosity. 

The interactions of the photons can be modelled by assuming that each photon can either 
interact directly or appear resolved through its fluctuations into hadronic states. In leading 
order QCD this model leads to three different event classes for 77 interactions: direct, single- 
resolved and double-resolved, where resolved means that the partons (quarks or gluons) inside 
the hadronic photon take part in the hard interaction. The probability to find partons in the 
photon is parametrised by parton density functions. 

At LEP the photons are radiated by the beam electrons^) and carry mostly small negative 
four- momenta squared, Q 2 . In this paper events are considered only if the electrons are scattered 
at small angles and are not detected. Both photons are therefore quasi- real (Q 2 « GeV 2 ). 
Differential hadron production cross-sections are measured as a function of the transverse mo- 
mentum and the pseudorapidity of charged hadrons. The inclusive cross-sections are compared 
to Monte Carlo (MC) event generators in leading order a s and theoretical calculations in next- 
to-leading order (NLO) a s for this process. 

1 In the OPAL coordinate system the x axis points towards the centre of the LEP ring, the y axis points 
upwards and the z axis points in the direction of the electron beam. The polar angle 9, the azimuthal angle <p 
and the radius r denote the usual spherical coordinates. The pseudorapidity rj is defined as 77 = — lntan(0/2). 

2 Positrons are also referred to as electrons. 
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2 The OPAL detector 



A detailed description of the OPAL detector can be found elsewhere [5]. The central tracking 
was performed inside a solenoidal magnet which provided a uniform axial magnetic field of 
0.435 T along the beam axis. Starting with the innermost components, the tracking system 
consisted of a high precision silicon microvertex detector, a precision vertex drift chamber, 
a large volume jet chamber with 159 layers of axial anode wires and a set of z chambers 
measuring the track coordinates along the beam direction. The transverse momenta, pr, of 
tracks are measured with a precision parametrised by <7 Pt /pt = y^.02 2 + (0.0015 • pr) 2 (pt 
in GeV) in the central region |cos#| < 0.73. In this paper transverse is always defined with 
respect to the beam axis. 

The magnet was surrounded in the barrel region (| cos6>| < 0.82) by a lead glass electro- 
magnetic calorimeter (ECAL) and a hadronic sampling calorimeter (HCAL), which in turn 
were surrounded by muon chambers. Similar layers of detectors were installed in the endcaps 
(0.82 < | cos#| < 0.98). The small angle region from 47 to 140 mrad around the beam pipe 
on both sides of the interaction point was covered by the forward calorimeters (FD) and the 
region from 33 to 59 mrad by the silicon tungsten luminometers (SW). The latter were used to 
determine the luminosity by counting small-angle Bhabha scattering events. 



3 Kinematics and Monte Carlo simulation 



The properties of the two interacting photons {i 
momentum transfers Q 2 and their invariant mass, 
angle Q\ relative to the beam direction by 



= 1,2) are described by their negative four- 
Each Q 2 is related to the electron scattering 



Qt = -(ft - V'if « 2£^(l - cos 00, (1) 

where pi and p\ are the four-momenta of the beam electron and the scattered electron, respec- 
tively, and Ei and E[ are their energies. Events with detected scattered electrons (single-tagged 
or double-tagged events) are excluded from the analysis. Driven by the angular acceptance of 
the FD and SW calorimeters a value of Q 2 = 4.5 GeV 2 is used in this analysis to separate the 
quasi-real photons of untagged events from the virtual photons of tagged events. The median 
Q 2 resulting from this definition cannot be determined from data since the scattered electrons 
are not tagged. For the kinematic range of this analysis the MC simulations predict the median 
Q 2 to be of the order of 10~ 4 GeV 2 . 

The hadronic invariant mass of the photon-photon system, W, can be obtained from the 
energies and momenta (£h, pt) of final state particles: 

^ 2 = (E^) 2 -(E^) 2 - (2) 

h h 

The MC generators PYTHIA 5.722 [6] and PHOJET 1.10 [7] have been used to simulate 
photon-photon interactions. PYTHIA uses the SaS-lD parametrisation [8] for the photon 
parton densities and PHOJET uses the GRV parametrisation [9]. 



5 



All relevant background processes were studied using MC generators. Multihadronic events 
(e + e~ — > qq(7)) were simulated with PYTHIA 6.125. KORALZ 4.02 [10] was used to generate 
the process e + e~ — > t + t~(7) and VERMASEREN [11] to generate e + e~ — > e + e - r + T~. Deep- 
inelastic e7 events were simulated with HERWIG 5.9 [12]. 

All signal and background MC samples were generated with full simulation of the OPAL 
detector [13]. They are analysed using the same reconstruction algorithms as for the data. 



4 Event selection 

The production of charged hadrons was studied using data taken at e + e~ centre-of-mass en- 
ergies, y/s^, from 183 to 209 GeV, amounting to a total integrated luminosity of 612.8 pb _1 . 
The luminosity weighted average centre-of-mass energy is 195.8 GeV. Two-photon events are 
selected with the following set of cuts: 

• At least six tracks must have been found in the tracking chambers. A track is required 
to have a minimum transverse momentum of 120 MeV with respect to the z axis and at 
least 40 hits in the central jet chamber. The number of measured hits in the jet chamber 
must be more than half of the number of possible hits given the track direction. The 
radial distance of nearest approach of the track to the primary vertex has to be less than 
0.15 cm. 

• The visible invariant hadronic mass calculated from the position and the energy of the 
clusters measured in the ECAL has to be greater than 3 GeV. 

• The sum of all energy deposits in the ECAL and the HCAL has to be less than 50 GeV 
to remove background from hadronic Z decays in events with a radiative return to the Z 
peak. 

• The missing transverse momentum of the event measured in the ECAL and the FD has 
to be less than 8 GeV. 

• To reject events with scattered electrons in the FD or SW, the total energy sum measured 
in the FD and SW has to be less than 60 GeV. 

• The background due to beam-gas or beam-wall interactions is reduced by requiring the 
radial distance of the primary vertex from the beam axis to be less than 2 cm and the 
distance from the nominal vertex position along the z direction to be less than 3 cm. 

After all cuts the data sample contains 1 144 035 events. The main remaining background 
processes are multihadronic Z decays, e + e~ — ► e + e~r + r~ and deep-inelastic e7 scattering. 
Other background processes are found to be negligible. The multihadronic background is 
mainly reduced by the cut on the sum of the energy measured by the ECAL and the HCAL 
and by the cut on the missing transverse momentum. The e + e~ — > e + e _ r + T~ processes are 
reduced by the cuts on the number of tracks and deep-inelastic e7 events are rejected by the 
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cut on the energy in SW and FD. The trigger efficiency for events in this region of phase space 
is close to 100% [4] and no correction is applied. 

From the MC simulations it is estimated that after all cuts the total remaining background 
is below 2% overall, but increasing to up to 50% at very high transverse momenta of the 
charged particles. The background is subtracted bin-by-bin from the distributions measured 
before corrections for detector resolution and acceptance are applied. The signal MC genera- 
tors PHOJET and PYTHIA are found to underestimate the cross-section. In particular they 
show a shape different from that in the data of the pp-distribution of charged hadrons. To 
study a potential bias resulting from this deficiency the MC events have been reweighted by a 
suitable p T -dependent function to resemble the data. Both the original and the reweighted MC 
distributions are used in the analysis and any differences seen are included in the systematic 
uncertainty. 



5 Analysis 

Only particles with a lifetime r > 0.3 ns are used to define the primary charged hadronic 
multiplicity. The primary charged hadrons originate either directly from the primary interaction 
or from the decay of particles with a lifetime r < 0.3 ns. The decay products of A and Kg 
particles are hence considered as primary hadrons in this analysis. The measured transverse 
momentum and pseudorapidity distributions of charged hadrons have to be corrected for the 
losses due to event and track selection cuts, the acceptance and the resolution of the detector. 
This is done using the MC events which were processed by the full detector simulation and 
reconstruction chain. The data are corrected by multiplying the experimental distribution, 
e.g. of the transverse momentum px, with correction factors which are calculated as the bin- 
by-bin ratio of the generated and the reconstructed MC distributions: 



da \ I dj>T / generated 

/ d(T 



MC 

(3) 



V / corrected ( ^ a ~\ V / measure d 

V dPT / reconstructed 

The ratio is calculated using both PYTHIA and PHOJET, with the mean value used to cor- 
rect the data. The pseudorapidity distributions are corrected in the same way. This method 
only yields reliable results if migration effects between bins due to the finite resolution of the 
measurement are small. The bins of the px and \rj\ distributions have therefore been chosen to 
be significantly larger than the resolution expected from the MC simulation. To avoid regions 
where the detector has little or no acceptance, all measurements of charged hadron distributions 
were restricted to the range \r]\ < 1.5. The average transverse momentum (px) and the average 
pseudorapidity (\r)\) in each bin is calculated directly from the data since detector corrections 
are small compared to the statistical uncertainties. 

In order to be able to use Eq. [3] for the detector correction in each bin of W, the bins in W 
must be larger than the experimental resolution and the average reconstructed hadronic invari- 
ant mass, (W rec ), should be approximately equal to the average generated hadronic invariant 
mass, (Wg en ). In this case bin-to-bin migrations are minimised. The visible invariant mass, 
H^vis, is determined from all tracks and calorimeter clusters, including the forward detectors and 
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the silicon tungsten calorimeters, after applying a matching algorithm to avoid double counting 
of particle momenta [15]. 

The ratio of W gen and Wis depends on the event kinematics, and was therefore determined 
separately in two distinct regions of phase space. The discriminating variable chosen for this 
purpose is pT,max, the highest transverse momentum of any track in the event. The average 
visible hadronic invariant mass, (Wvis), and the resolution on Wd S as a function of the generated 
hadronic invariant mass W gen for events with px.max smaller and larger than 5 GeV are shown 
in Fig. [0(a) and (c). The polynomial fits superimposed on Fig. [TJb) and (d) are used as a 
correction function so that (W gen )/ (Wee) ~ 1- 

The distributions are measured for invariant masses 10 < W < 30 GeV, 30 < W < 50 GeV, 
50 < W < 125 GeV and 10 < W < 125 GeV where W represents the hadronic invariant mass 
after correcting for detector effects. To facilitate comparisons to the results in [1], the ranges 
W > 30 GeV and W > 50 GeV are also considered. 



6 Systematic uncertainties 

The following sources of systematic uncertainties have been investigated: 

• Significant amounts of background need to be subtracted from the data at large pp, 
where the dominant source of background is multihadronic events (e + e~ — > qq(7))- The 
amount of this background subtracted from the data is varied by ±10% to estimate the 
uncertainty associated to this procedure. The amount by which the background must 
be varied is estimated from studying the MC description of the data in regions of the 
phase space where the background dominates. The effect on the measured cross-section 
is usually 1 % or less for low transverse momenta, increasing to up to 7 % at the highest 
momenta measured. 

• The selection criteria described in Section H] are varied simultaneously both to be more 
restrictive and to allow more events into the analysis to exclude a strong dependence on 
the event selection. Selection criteria based on energy measurements are varied by 5% 
in the ECAL and HCAL, and by 10% in the FD and SW calorimeters. The number of 
tracks required is changed by ±1. The allowed radial distance of the tracks is varied by 
5%. The uncertainty on the cross-section derived from all these variations is typically 
1-6%. 

• Both PHOJET and PYTHIA have been used to calculate the correction factors applied 
to the data. The resulting distributions are averaged and the difference between the two 
distributions is used to define the systematic uncertainty. The distributions obtained from 
both programs have been reweighted for a better description of the data. The difference 
between using the reweighted and the original MC distributions to calculate the correction 
factors has been included in the systematic uncertainty. The uncertainty derived from 
this study is below 5%. 
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• The systematic uncertainty due to the energy scale of the ECAL for the range of energies 
in this analysis was estimated by varying the reconstructed ECAL energy in the MC by 
±3% [14]. The cross-sections change by up to 4% due to this variation. 

• Studying vertex and net charge distributions it is estimated that about 2% of the selected 
events are due to beam-gas or beam-wall interactions. This remaining background is 
treated systematic uncertainty of 2%. 

Systematic uncertainties due to the modelling of the detector resolution for the measurement 
of tracks [4] and due to the luminosity measurement were found to be negligible. The total 
systematic uncertainty was obtained by adding all systematic uncertainties in quadrature. 



7 Results and Conclusions 

The differential inclusive cross-sections dcr/dpr for charged hadrons in four ranges of the 
hadronic invariant mass W are shown in Fig. [2] and are given in Table [1] and Table [2j The 
data points in the figure indicate the bin centre. Both data and calculations are presented for 
quasi-real photons of virtualities Q 2 < 4.5 GeV 2 , as described in Sect. El 

Calculations in NLO QCD [16] are compared to the data. The cross-sections are calculated 
using the QCD partonic cross-sections in NLO for direct, single- and double-resolved processes. 
The hadronic cross-section is a convolution of the Weizsacker- Williams effective photon distri- 
bution, the parton distribution functions and the fragmentation functions of [17]. The AFG-HO 
parametrisation of the parton densities of the photon [18] is used with =221 MeV. The 
renormalisation and factorisation scales in the calculation are set equal to pp- The cross-section 
calculation was repeated for the kinematic conditions of the present analysis. For the differ- 
ential cross-section dcr/dpr a minimum pr of 1.5 GeV is required to ensure the validity of the 
perturbative QCD calculation. Even at pr = 1-5 GeV the cross-sections change by up to 80% 
when varying the renormalisation and factorisation scales by factors of two. This uncertainty 
decreases rapidly to between 10% and 15% for p T = 3.5 GeV and above. The differential 
cross-section do"/d|?]| is hence restricted to the region p^ > 3.5 GeV to allow a meaningful 
comparison. 

In Fig. [2] the NLO calculation lies significantly below the data for transverse momenta 
greater than about 10 GeV, which can be reached in the highest W range only. The predictions 
of PHOJET and PYTHIA at high px (not shown) are similar to the NLO calculation. 

The differential cross-section is nearly independent of \rj\ in the range measured as can be 
seen in Fig. [3j The data points in the figure indicate the bin centre. The cross-section values 
are given in Table [31 The NLO calculation reproduces the data well within the uncertainties of 
the calculation. PYTHIA and PHOJET describe the shape of the distributions correctly for all 
regions of the phase space measured, but are below the data in normalisation for large values 
of W. 

Fig. H](a) and (b) and TableHlshow the differential cross-section dcx/dpr for charged hadrons 
for W > 30 GeV and W > 50 GeV to facilitate a comparison with a recent measurement by L3 
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of charged pions in the pseudorapidity range \rj\ < 1.0 [1]. The OPAL data shown in Fig. H] (c) 
and (d) and Table have been scaled to account for the reduced \rj\ range and for the fraction 
of charged pions of all charged hadrons using MC simulations. While the OPAL data points 
in Fig. H] (a) and (b) are plotted at the bin centre, in Fig. [4] (c) and (d) they are plotted at 
the mean transverse momentum across the bin, as for the L3 data. From this comparison it 
is evident that the distributions measured by OPAL fall more rapidly towards high transverse 
momenta than those measured by L3, leading to a disagreement between the two experiments 
at high transverse momenta and a better description of the OPAL data by NLO QCD than is 
the case for the L3 data. 
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-1.80 


1.69 


(1.92 ±0.04 ±0.13} 
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1.69 


(9.13±0.21±0.32} 


xlO 1 


1.80- 


-2.00 


1.89 


(1.13±0.03±0.08} 


xlO 2 


1.89 


(5.32±0.15±0.19} 


xlO 1 


2.00- 


-2.20 


2.09 


(6.68 ±0.23 ±0.44} 


xlO 1 


2.09 


(3.34 ±0.12 ±0.11} 


xlO 1 
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-2.40 


2.29 


(4.17±0.19±0.26} 


xlO 1 


2.29 


(2.22 ±0.10 ±0.07} 


xlO 1 


2.40- 


-2.60 


2.50 
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xlO 1 


2.50 


(1.47 ±0.08 ±0.04} 


xlO 1 


2.60- 


-2.80 


2.70 


(1.98±0.13±0.12} 


xlO 1 


2.70 


(1.00±0.06±0.03} 


xlO 1 
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-3.00 


2.90 


(1.45 ±0.12 ±0.09} 


xlO 1 


2.90 


(7.48 ±0.54 ±0.26} 
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xlO" 1 


5.43 
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xlO" 1 
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-8.00 


6.74 ±0.01 


(1.37±0.31±0.22} 


xlO" 1 
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xlO" 1 
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9.52 ±0.12 


(1.90 ±0.86 ±0.33} 


xlO' 2 


9.62 ±0.03 


(2.23 ±0.60 ±0.36} 


xlO" 2 



Table 1: Differential inclusive charged hadron production cross-sections dcr/dpr for < 1.5 
and in the W ranges 10 < W < 30 GeV and 30 < W < 50 GeV. The first uncertainty is the 
statistical uncertainty and the second uncertainty is the systematic uncertainty. No value is 
given if the error on (px) is less than 0.01. 
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50 < W < 125 GeV 


10 < W < 125 GeV 


Pt 


(Pt) 


dcr/dpr 




(Pt) 




da/dpT [pb/GeV] 


[GeV] 


[GeV] 


[pb/GeV] 




[GeV] 




[pb/GeV] 




0.12- 


-0.28 


0.20 


(1.00 ±0.00 ±0.03) 


xlO 4 


0.20 




(5.93±0.00±0.21) 


xlO 4 


0.28- 


-0.44 


0.35 


(7.77±0.02±0.22) 


xlO 3 


0.35 




(4.68 ±0.00 ±0.18) 


xlO 4 


0.44- 


-0.60 


0.51 


(4.58 ±0.01 ±0.13) 


xlO 3 


0.51 




(2.71±0.00±0.10) 


xlO 4 


0.60- 


-0.80 


0.69 


(2.37±0.01±0.07) 


xlO 3 


0.69 




(1.34 ±0.00 ±0.05) 


xlO 4 


0.80- 


-1.00 


0.89 


(1.12±0.01±0.03) 


xlO 3 


0.89 




(5.87±0.01±0.25) 


xlO 3 


1.00- 


-1.20 


1.09 


(5.56 ±0.04 ±0.14) 


xlO 2 


1.09 




(2.68 ±0.01 ±0.14) 


xlO 3 


1.20- 


-1.40 


1.29 


(2.92 ±0.03 ±0.08) 


xlO 2 


1.29 




(1.29 ±0.01 ±0.07) 


xlO 3 


1.40- 


-1.60 


1.49 


(1.66±0.02±0.05) 


xlO 2 


1.49 




(6.78 ±0.04 ±0.38) 


xlO 2 


1.60- 


-1.80 


1.69 


(9.74 ±0.15 ±0.30) 


xlO 1 


1.69 




(3.78 ±0.03 ±0.20) 


xlO 2 


1.80- 


-2.00 


1.89 


(6.10±0.11±0.21) 


xlO 1 


1.89 




(2.27±0.02±0.12) 


xlO 2 


2.00- 


-2.20 


2.09 


(4.01 ±0.09 ±0.15) 


xlO 1 


2.09 




(1.41 ±0.02 ±0.07) 


xlO 2 


2.20- 


-2.40 


2.29 


(2.53±0.07±0.10) 


xlO 1 


2.29 




(9.03 ±0.16 ±0.44) 


xlO 1 


2.40- 


-2.60 


2.50 


(1.76 ±0.06 ±0.09) 


xlO 1 


2.50 




(6.18±0.13±0.31) 


xlO 1 


2.60- 


-2.80 


2.70 


(1.22±0.05±0.06) 


xlO 1 


2.70 




(4.26 ±0.11 ±0.22) 


xlO 1 


2.80- 


-3.00 


2.90 


(9.52 ±0.44 ±0.50) 


xlO 


2.90 




(3.25±0.09±0.17) 


xlO 1 


3.00- 


-3.50 


3.23 


(5.99 ±0.22 ±0.30) 


xlO 


3.23 




(1.93 ±0.05 ±0.10) 


xlO 1 


3.50- 


-4.00 


3.73 


(3.33±0.16±0.22) 


xlO 


3.73 




(1.03±0.03±0.06) 


xlO 1 


4.00- 


-5.00 


4.40 


(1.52±0.08±0.15) 


xlO 


4.43 




(4.66 ±0.16 ±0.34) 


xlO 


5.00- 


-6.00 


5.43 


(7.02 ±0.83 ±0.81) 


xlO" 1 


5.43 




(1.71±0.11±0.13) 


xlO 


6.00- 


-8.00 


6.83 


(3.89 ±0.45 ±0.49) 


xlO" 1 


6.79 




(7.28 ±0.53 ±0.62) 


xlO" 1 


8.00- 


15.00 


10.18±0.01 


(8.40 ±1.34 ±0.90) 


xlO^ 2 


10.00 




(1.14±0.14±0.09) 


xlO" 1 


15.00- 


-25.00 


18.26 ±0.09 


(2.46 ±0.96 ±0.26) 


xlO" 2 


18.27±0.08 


(2.76 ±0.99 ±0.22) 


xlO" 2 



Table 2: Differential inclusive charged hadron production cross-sections dcr/dpr for M < 1-5 
and in the W ranges 50 < W < 125 GeV and 10 < W < 125 GeV. The first uncertainty is 
the statistical uncertainty and the second uncertainty is the systematic uncertainty. No value 
is given if the error on (px) is less than 0.01. 
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10 < W < 30 GeV 


30 < W < 50 GeV 


\v\ 


<M> 


der/d|?y| [pb] 


(N) 


dcr/d?? [pb] 


0.0-0.3 


0.149 


3.26±0.45±0.29 


0.151 


2.49±0.26±0.07 


0.3-0.6 


0.445 


3.52±0.50±0.39 


0.451 


2.41±0.25±0.08 


0.6-0.9 


0.754 


3.11±0.47±0.41 


0.759 


2.48±0.26±0.07 


0.9-1.2 


1.047 


2.58±0.43±0.24 


1.049 


2.15±0.26±0.17 


1.2-1.5 


1.348 


2.40±0.44±0.28 


1.341 


1.74±0.23±0.07 





50 < W < 125 GeV 


10 < W < 125 GeV 




7/ 




( 


T) 


> 


da/dM [pb] 


( 


7/ 


) 


dcr/d?, [pb] 


0.0-0.3 
0.3-0.6 
0.6-0.9 
0.9-1.2 
1.2-1.5 


0.149 
0.449 
0.759 
1.048 
1.345 


4.17±0.27±0.63 
3.63±0.25±0.53 
3.57±0.25±0.46 
2.81±0.23±0.34 
2.81±0.22±0.25 


0.149 
0.448 
0.757 
1.047 
1.345 


10.90±0.65±0.90 
9.92 ±0.62±0.65 
9.89 ±0.63±0.63 
8.04 ±0.58±0.67 
7.42 ±0.55±0.44 



Table 3: Differential inclusive charged hadron production cross-sections d<r/d|?7| for 
p T > 3.5 GeV and in the W ranges 10 < W < 30 GeV, 30 < W < 50 GeV, 50 < W < 125 GeV 
and 10 < W < 125 GeV. The first uncertainty is the statistical uncertainty and the second 
uncertainty is the systematic uncertainty. 
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W > 30 GeV 


W > 50 GeV 


Pt 


(Pt) 


dcr/dpr 




(Pt) 


da/dp T 




[GeV] 


[GeV] 


[pb/GeV] 




[GeV] 


[pb/GeV] 




0.12- 


-0.28 


0.20 


(2.25 ±0.00 ±0.06} 


xlO 4 


0.20 


(1.10±0.00±0.04) 


xlO 4 


0.28- 


-0.44 


0.35 


(1.75 ±0.00 ±0.05) 


xlO 4 


0.35 


(8.50 ±0.02 ±0.28) 


xlO 3 


0.44- 


-0.60 


0.51 


(1.03 ±0.00 ±0.03) 


xlO 4 


0.51 


(5.02±0.02±0.16) 


xlO 3 


0.60- 


-0.80 


0.69 


(5.31 ±0.02 ±0.13) 


xlO 3 


0.69 


(2.61±0.01±0.08) 


xlO 3 


0.80- 


-1.00 


0.89 


(2.46 ±0.01 ±0.06) 


xlO 3 


0.89 


(1.24 ±0.01 ±0.04) 


xlO 3 


1.00- 


-1.20 


1.09 


(1.20±0.01±0.03) 


xlO 3 


1.09 


(6.17 ±0.05 ±0.18) 


xlO 2 


1.20- 


-1.40 


1.29 


(6.19 ±0.05 ±0.16) 


xlO 2 


1.29 


(3.25 ±0.04 ±0.09) 


xlO 2 


1.40- 


-1.60 


1.49 


(3.44 ±0.04 ±0.09) 


xlO 2 


1.49 


(1.85±0.03±0.05) 


xlO 2 


1.60- 


-1.80 


1.69 


(2.00 ±0.03 ±0.06) 


xlO 2 


1.69 


(1.09 ±0.02 ±0.03) 


xlO 2 


1.80- 


-2.00 


1.89 


(1.22 ±0.02 ±0.04) 


xlO 2 


1.89 


(6.87±0.18±0.20) 


xlO 1 


2.00- 


-2.20 


2.09 


(7.86 ±0.19 ±0.24) 


xlO 1 


2.09 


(4.51 ±0.14 ±0.14) 


xlO 1 


2.20- 


-2.40 


2.29 


(5.09±0.15±0.16) 


xlO 1 


2.29 


(2.86±0.11±0.10) 


xlO 1 


2.40- 


-2.60 


2.50 


(3.48 ±0.12 ±0.13) 


xlO 1 


2.50 


(2.01 ±0.09 ±0.08) 


xlO 1 


2.60- 


-2.80 


2.70 


(2.38 ±0.10 ±0.09) 


xlO 1 


2.70 


(1.38±0.08±0.06) 


xlO 1 


2.80- 


-3.00 


2.90 


(1.83 ±0.09 ±0.07) 


xlO 1 


2.90 


(1.08±0.07±0.05) 


xlO 1 


3.00- 


-3.50 


3.23 


(1.12 ±0.04 ±0.04) 


xlO 1 


3.23 


(6.78 ±0.34 ±0.27) 


xlO 


3.50- 


-4.00 


3.73 


(6.08 ±0.32 ±0.25) 


xlO 


3.73 


(3.72 ±0.26 ±0.17) 


xlO 


4.00- 


-5.00 


4.43 


(2.89±0.16±0.16) 


xlO 


4.43 


(1.74 ±0.12 ±0.13) 


xlO 


5.00- 


-6.00 


5.43 


(1.35 ±0.14 ±0.08) 


xlO 


5.43 


(8.23 ±1.22 ±0.73) 


xlO" 1 


6.00- 


-8.00 


6.80 


(6.61 ±0.79 ±0.48) 


xlO" 1 


6.83 


(4.57 ±0.72 ±0.44) 


xlO" 1 


8.00- 


15.00 


10.11 


(1.22 ±0.20 ±0.09) 


xlO" 1 


10.25 


(1.04 ±0.21 ±0.10) 


xlO" 1 


15.00- 


-25.00 


18.16±0.03 


(3.43 ±2.46 ±0.27) 


xlO" 2 


18.14±0.03 


(3.43 ±1.34 ±0.32) 


xlO" 2 



Table 4: Differential inclusive charged hadron production cross-sections da/dp^ for < 1.5 
and in the ranges W > 30 GeV and W > 50 GeV. The first uncertainty is the statistical 
uncertainty and the second uncertainty is the systematic uncertainty. No value is given if the 
error on (p T ) is less than 0.01. 
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W > 30 GeV 


W > 50 GeV 


Pt 


(Pt) 


dcr/dpr 




(Pt) 


dcr/dpr 




[GeV] 


[GeV] 


[pb/GeV] 




[GeV] 


[pb/GeV] 




0.12- 


-0.28 


0.20 


(1.38 ±0.00 ±0.04] 


xlO 4 


0.20 


(6.71 ±0.02 ±0.22) 


xlO 3 


0.28- 


-0.44 


0.35 


(1.05±0.00±0.03) 


xlO 4 


0.35 


(5.11 ±0.01 ±0.17) 


xlO 3 


0.44- 


-0.60 


0.51 


(5.85 ±0.01 ±0.15) 


xlO 3 


0.51 


(2.85 ±0.01 ±0.09) 


xlO 3 


0.60- 


-0.80 


0.69 


(2.80 ±0.01 ±0.07) 


xlO 3 


0.69 


(1.37 ±0.01 ±0.04) 


xlO 3 


0.80- 


-1.00 


0.89 


(1.20±0.01±0.03) 


xlO 3 


0.89 


(6.04 ±0.04 ±0.19) 


xlO 2 


1.00- 


-1.20 


1.09 


(5.55 ±0.04 ±0.14) 


xlO 2 


1.09 


(2.83 ±0.03 ±0.08) 


xlO 2 


1.20- 


-1.40 


1.29 


(2.75 ±0.03 ±0.07) 


xlO 2 


1.29 


(1.43 ±0.02 ±0.04) 


xlO 2 


1.40- 


-1.60 


1.49 


(1.49 ±0.02 ±0.04) 


xlO 2 


1.49 


(7.88 ±0.14 ±0.20) 


xlO 1 


1.60- 


-1.80 


1.69 


(8.51 ±0.14 ±0.24) 


xlO 1 


1.69 


(4.58±0.10±0.12) 


xlO 1 


1.80- 


-2.00 


1.89 


(5.04±0.11±0.15) 


xlO 1 


1.89 


(2.77 ±0.08 ±0.08) 


xlO 1 


2.00- 


-2.20 


2.09 


(3.25±0.09±0.10) 


xlO 1 


2.09 


(1.80±0.06±0.06) 


xlO 1 


2.20- 


-2.40 


2.29 


(2.10 ±0.07 ±0.07) 


xlO 1 


2.29 


(1.15 ±0.05 ±0.04) 


xlO 1 


2.40- 


-2.60 


2.50 


(1.41 ±0.05 ±0.05) 


xlO 1 


2.50 


(7.90 ±0.41 ±0.33) 


xlO 


2.60- 


-2.80 


2.70 


(9.75 ±0.45 ±0.39) 


xlO 


2.70 


(5.43 ±0.33 ±0.23) 


xlO 


2.80- 


-3.00 


2.90 


(7.49 ±0.40 ±0.29) 


xlO 


2.90 


(4.19 ±0.30 ±0.18) 


xlO 


3.00- 


-3.50 


3.23 


(4.72 ±0.20 ±0.18) 


xlO 


3.23 


(2.79±0.16±0.11) 


xlO 


3.50- 


-4.00 


3.73 


(2.54±0.15±0.10) 


xlO 


3.73 


(1.53 ±0.12 ±0.07) 


xlO 


4.00- 


-5.00 


4.43 


(1.19 ±0.07 ±0.06) 


xlO 


4.43 


(7.03±0.56±0.51) 


xlO" 1 


5.00- 


-6.00 


5.43 


(5.99 ±0.68 ±0.38) 


xlO" 1 


5.43 


(3.66 ±0.58 ±0.32) 


xlO" 1 


6.00- 


-8.00 


6.80 


(2.85 ±0.36 ±0.21) 


xlO" 1 


6.83 


(1.99 ±0.33 ±0.19) 


xlO" 1 


8.00- 


15.00 


10.11 


(5.27±0.93±0.40) 


xlO" 2 


10.25 


(4.28 ±0.91 ±0.40) 


xlO" 2 


15.00- 


-25.00 


18.16±0.03 


(1.78 ±1.32 ±0.14) 


xlO" 2 


18.14±0.03 


(1.80 ±0.78 ±0.17) 


xlO" 2 



Table 5: Differential inclusive charged pion production cross-sections da/dpr for \r]\ < 1.0 and 
in the ranges W > 30 GeV and W > 50 GeV. The first uncertainty is the statistical uncertainty 
and the second uncertainty is the systematic uncertainty. No value is given if the error on (pr) 
is less than 0.01. 
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Figure 1: The correlation between the generated hadronic invariant mass Wgen and the visible 
hadronic invariant mass W v i s in two regions of pr,max, the maximum transverse momentum 
of any track in the event, for PHOJET and PYTHIA MC events. The vertical bars show 
the standard deviation in each bin in (a,c) and the uncertainty on the mean in (b,d) where 
larger than the marker size. The horizontal lines indicate the bin width. The polynomial fit 
determines the correction function for W v ; s . 
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Figure 2: Differential inclusive charged hadron production cross-sections, da/dp^, for \r]\ < 1.5 
in the W ranges (a) 10 < W < 30 GeV; (b) 30 < W < 50 GeV; (c) 50 < W < 125 GeV 
and (d) 10 < W < 125 GeV. The error bars show the statistical and systematic uncertainties 
added in quadrature when larger than the marker, which indicates the bin centre. The data 
are compared to an NLO calculation [16]. 
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Figure 3: Differential inclusive charged hadron production cross-sections, da/d^l, for 
p T > 3.5 GeV and in the W ranges (a) 10 < W < 30 GeV; (b) 30 < W < 50 GeV; (c) 
50 < W < 125 GeV and (d) 10 < W < 125 GeV. The data are compared to the PHOJET 
and PYTHIA models. The data are also compared to an NLO calculation [16]. The error bars 
show the statistical and systematic uncertainties added in quadrature. The markers are shown 
at the centre of the bin. 
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Figure 4: Differential inclusive charged hadron production cross-sections, dcr/dpx, f° r \v\ < 1-5 
and in the W ranges (a) W > 30 GeV and (b) W > 50 GeV. The differential inclusive charged 
pion production cross-sections da/dp T for \rj\ < 1.0 in the W ranges (c) W > 30 GeV and 
(d) W > 50 GeV are inferred from the charged hadron cross-sections using MC information to 
facilitate the comparison to the L3 measurements [1]. The error bars show the statistical and 
systematic uncertainties added in quadrature when larger than the marker. While the OPAL 
data points in (a) and (b) are plotted at the bin centre, in (c) and (d) they are plotted at the 
mean transverse momentum across the bin, as for the L3 data. 
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